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Abstract 

It is siiown tliat tlie gluons produced witli small transverse mo- 
menta k±_ inside a jet {k± < A;^*, /c^"* ~^ 0) ^ire independently emitted 
from the primary parton, as QCD coherence suppresses their show- 
ering. Consequently, the low k±_ gluons follow a Poisson distribution, 
very much like the soft photons radiated by a charged particle in QED. 
On the contrary, the distribution of gluons with limited absolute mo- 
menta \k\ < k^'^^ remains non- Poissonian even for small A;'^"*. It will 
be interesting to find out to what extent this perturbative prediction 
for partons survives the hadronization process. 
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1 Introduction 



It is by now well established that the multiplicity distributions of hadrons 
produced in high energy collisions are substantially wider than the Poisson 
distribution. This can be understood in a large class of models based on 
branching processes |]T|, the quark gluon cascade in perturbative QCD be- 
ing a relevant example @J^. In this paper we point out that the specific 
properties of the QCD cascade allow one to select simple regions of phase 
space where the produced partons are less correlated and show a narrow 
multiplicity distribution. The experimental check of our predictions would 
therefore provide a novel test of the role of perturbative QCD in multiparticle 
production. 

It is convenient to characterize the multiplicity distribution of particles 
in an event by its global factorial moments f^'^^ =< n{n — l)...{n — q + l) > 
of order g, or, written in terms of the g-particle inclusive density 

f''\Q)= I p^''Hh...kg;Q)dh...dk„ (1) 

with an integral over the full phase space available at the energy scale Q for 
the process. The ratios F*^'^^ = f^^^ /n^ were extensively studied in perturba- 
tive QCD to various degrees of accuracy 0,^0. The asymptotic results are 
obtained in the Double Logarithmic Approximation (DLA) and the normal- 
ized moments approach a limit F^'^\Q) Fg, corresponding to KNO mul- 
tiplicity scaling |jl|J^. The numbers Fg are independent of the coupling, thus 
providing a universal characteristic of the multiplicity distribution. These 
moments essentially differ from = 1 as for the Poisson distribution and 
then the asymptotic KNO distribution essentially differs from the Dirac delta 
function. 

In analogy with (0) we define the cut moments 

Q)= [ p^'^\h ■■■K; Q)dh . . . dkg, (2) 

where the phase space integration is now restricted by Tc{C, Q) with a cut 
variable C . In this paper we discuss two types of cuts to be applied to all 
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particles in the final state: the momentum cut \ki\ < /c'^"*, and the trans- 
verse momentum cut k±^i < kf^^, for i = l,...,q. Clearly, the moments (^) 
determine the multiplicity distribution of particles produced in the restricted 
phase space, and as such they provide a more differential characteristics than 
the global quantities. For the maximal /c'^"* and /c^"* at a given energy scale 
Q the global quantities are retained. 

Depending on the choice of the cut, the moments in Eq.(0) are probing 
quite different physical phenomena. Of special importance in this discussion 
is the coherence of the soft gluon emission from the harder partons which 
yields an angular ordering of subsequent emissions along the branching pro- 
cess Q. In particular, this property leads to the suppression of the soft 
particle production and to a fiat, energy independent rapidity distribu- 
tion of particles with low transverse momenta in the parton jet, reminiscent 



of the QED Bremsstrahlung |T0 



To see the effect of the angular ordering on our cut moments consider 
first the momentum cut. Since the dominant source of particle production 
in the QCD cascade is the gluon splitting, the existence of one gluon simply 
enhances the probability to find another one nearby in momentum space 
and this yields positive correlations which do not vanish even for a small 
momentum cut k'^^^. Such correlations would also persist with the k±_ cutoff if 
there was no angular ordering. With the angular ordering taken into account 
however, a cut in the transverse momentum of one parton yields an even 
stronger restriction of the available phase space for the next parton in the 
cascade, and consequently, the correlations are strongly damped near the 
low k± threshold. Because of this suppression of secondary emissions in the 
k± cut phase space one obtains ultimately the independent emission of low 
k± gluons and therefore a Poisson distribution very much in analogy to the 
multiple soft photon Bremsstrahlung in QED. 

Differences between the momentum cut and k± cut moments become 
more pronounced towards small cutoffs. Fortunately, this is the region where 
the simple approximations are most reliable. We derive here the evolution 
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equations for both types of moments, which can be solved recursively to any 
order of the perturbation theory within the DLA. These results provide the 
qualitatively correct picture as can also be inferred from the Monte Carlo 
simulations at the parton level. 

It will be interesting to find out whether the experimental data indeed 
show the approach to a Poisson distribution as one might expect from the Lo- 
cal Parton Hadron Duality (LPHD) picture, originally proposed for single in- 
clusive distributions , but successful also in a broader range of phenomena 



TT|. Alternatively, the hadronization process could destroy the coherence 
properties of the parton cascade; this happens in a specific hadronization 
model studied here for small fc'^"*, /c^"* < 1 GeV. Hence the observation of the 
above regularities in the hadronic final state would provide us with a clue on 
the hadronization mechanism. 

After presenting a few basic relations of the DLA scheme in Chapter 2, we 
derive in Chapter 3 the threshold behaviour of both families of cut moments 
in the leading perturbative order. In Chapter 4 recursive evolution equations 
are derived which generate perturbative expansions to arbitrary order for 
constant and running Ug- In Chapter 5 we discuss for comparison the Monte 
Carlo results and conclude with Chapter 6. The Appendices contain details 
of various derivations and a sample of explicit expressions for cut moments 
up to the third non-leading correction. 



2 Particle correlations in jets in DLA 

In this paper we consider multiplicity distributions of hadrons or partons in 
a jet of primary energy P and half opening angle 6 in limited regions of 
phase space. The multiparton correlations in a jet evolving from the primary 
parton a {a = q, g) can be derived from the generating functional Zp^a which 
in DLA obeys the following integral equation |^ 

Zp^,{u} = exp f / Mp,a{K)[u{K)ZK,g{u} - l]d'K] . (3) 
\Jrp(K) J 
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This equation refers to the multiparton final state excluding the leading par- 
ticle. In the subscript P denotes collectively the momentum vector of 
the parent parton and the half opening angle (P = {P, 0}) of the jet it 
generates; Tp{K) stands for the phase space of the intermediate parent 
Vp{K) = {K : K < P, Qkp < ©, KQkp > Qo}, where Qo is a transverse mo- 
mentum cutoff parameter). The probability Ai p,a{K) for the bremsstrahlung 
of a single gluon off the primary parton a reads for small angle B pk 



Mp,a{K)d^K = Cail{KQpK) — — (4) 

where Cg = 1, = Cf/A^^c = 4/9 for initial gluons and quarks respectively. 
The multiplicity anomalous dimension is related to the strong coupling by 
7q = 2Ncas/'K and is taken as 7o(y) = P'^/iy + A) where y = ln(A;_|_/A), 
A = In ^ and = ANc / {^^c ~\^f)'i here Nc and nj denote the numbers 
of colours and flavours and A is the QCD scale. 

The inclusive densities are obtained from Taylor expansion of the gener- 
ating functional around u{k) = 1 

pp^aih, K) = S''Zp^a{u}/Su{ki)...6u{kg) (5) 

likewise the connected (cumulant) correlation functions from 

r^^i(fci,... ,kq) = 5'^\nZp,a{u}/Su{k,) . . .Su{kq) \u=i . (6) 

These differential distributions satisfy the following integral equations |14 

rS?_i(A;i, ...,k,)= S^l^^k,, ...,kq) + J d'KMp,a{K)p^;^l{K ...,kq), 

(7) 

dp!a,nestiku ■ ■ ■ , kg) = A^P.alA^Opi^/^ (^2, ■ ■ ■ , kg) + cycl. (8) 

The cut moments fj:'^^ are obtained as integrals over the particle density 
distributions p*^^-* as in eq. (0), likewise the cumulant moments c'-'^^ as integrals 
over the connected correlation functions T^'^K 

The global moments integrated over full phase space can be derived more 
simply from the integrated generating function Z{Y,u) which - in analogy 
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to eq. d^) - fulfils the evolution equation ||T3| 

Z,{Y,u) = exp{ r dyCa^l{y){Y - y)[uZ{y,u) - 1]} (9) 
Jo 

where the evolution variable in small angle approximation is y = hi^. 
The global moments are then derived from /^''^ = df' Z / du'^\u=i and c^''^ = 
d'^hiZ / du''\u=i and we recall that these moments are related by (see, for 
example |]T2[) 

j{q+l) = ( ^\f{g-m)^{m+l) /^gN 

Frequently used are the moments F^'^^ = f^'^^/n'^ and K^'^^ = c^'^'> /n'' normal- 
ized by the mean multiplicity n = /"^^^ = c'-^-*. 

For later use we also present here the results for the global factorial mo- 
ments /^'^"^^ to the lowest perturbative order, an approximation suitable for 
small Y. To this end, one obtains first, after q-fold differentiation of InZ 
over M in (^, the equation 

c^\Y) = r dyCall{y){Y -y){qf'^'^\y) + f'^\y)). (11) 



Then, starting with f^^^ = 1, one finds the higher moments by iteration 
whereby in leading order of Y the second term with f^'^^ in ([TT|) can be 
neglected against qf^"^^^- After inserting (0) into (|T0|) one obtains finally 
for small Y 

f^^\Y)^f,cl^l'^Y'^, (12) 

thus the moments vanish at the threshold Y ^ [PQ —>■ Qq). The coeffi- 
cients fg in (p!2| ) are given recursively by 

f - fkfq-k-l ^\ f^o\ 

and the first coefficients read 

1 1 17 

/o = 1, /i = 2' /a = 2' /s = ^5 • • • ■ (14) 

Note that the relation (0) is different from the analogous relation which 
determines global moments at large Y ( P^JI^ ). 
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3 Multiplicity Moments for small cut-off 



In this section we derive the multiphcity moments for the phase space with 
cyhndrical cut k± < A;™* and with spherical cut k < /c*^"* for small parameters 
/cj"* and /c*^"*. In these limiting cases various simplifications can be apphed: 
first, it is enough to take into account the terms of lowest order in the cou- 
pling. We show in the next section that terms of higher order are indeed 
suppressed for small transverse momenta. Secondly, in this limit, the solu- 
tions for running coupling approach those for fixed coupling, in which case 
close expressions can be obtained. 

3.1 Cylindrical phase space 

For the multiphcity in the /cj"* cyhnder (not including the leading particle of 
momentum P) we obtain 

n{kr, P, e) = ^I'oik^) I dy^l{y){Y - y), 

JQo k± JkjL P W JQ 

(15) 

where X±_ = In Now it can easily be seen that for X± <^ A the variation 
of the coupling can be neglected and ■Jq = For fixed coupling we find 

in lowest order 

n{X^,Y)=^',X4Y-Xj2), (16) 

which reduces for X± <^ y to n{X±,Y) — ^qX±Y corresponding to a con- 
stant particle density in \nk± and In A;. 

Here and in the rest of this section we write for simplicity the results 
for gluon jets. For quark jets the coupling 7q should be replaced by Cq^jg 
everywhere as we deal only with gluon emission from the primary parton. 

The higher order moments are obtained from (7), where we neglect the 
second term which is of higher order in as- The integration of F^^^ and 
over the momenta yields the cumulant and factorial moments in the 
given region of phase space. For ki kept fixed the integral over k2 - ■ ■ kg in 
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(8) yields the global factorial moment /'^''^ at scale k±i = ki'di. because of 
angular ordering the angles in all later emissions must obey '(92, ■ ' ' I'^q < "^i- 
So one finds 



c(^\X^,Y) = q f''"dyj',{y){Y-y)f(^-'\y). (17) 
Jo 

Hence in this approximation the cut moments are determined by the global 
ones. Again one recovers the fixed as limit for small Xj_. 

For the normalized cumulant moments in the case of fixed coupling we 
find from (0) and ( p!7D after inserting the moments in the approximation 



and for the case of small cut-off fc]"* we are interested in, i.e. for X± ^ Y 

ir(«)(X„F) = ^(^)"\ (19) 

Therefore, in this approximation the cumulants K^*^^ quickly decrease with 
the order q. It is interesting to note that this power like dependence on the 



order q corresponds to the "linked pair" ansatz for correlations |jT5[, usually 
written as K^'^'' = Ag[K^'^^]'^^^ , which is made in models where multi-particle 
correlations are built from 2-particle correlations. The higher cumulants in 



([T^) rise more slowly than in case of the negative binomial distribution in 



which case Ag = {q — 1)\ [jl6|. 



Finally, the factorial moments are obtained by solving (|ToD approximately 
through = 1 + g(g - l)ir(2) /2 + 0{K^^'^), which yields 

F(^)(Xx,r)-l + ^^^^. (20) 

o Y 

So we obtain the remarkable result that for small transverse momentum 
cut-off all factorial moments approach unity and therefore the multiplicity 
distribution becomes Poissonian. This is a consequence of the dominance 
of the single soft gluon emission at small k±, i.e., the absence of branching 
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processes with secondary gluon emissions. This behaviour is just analogous 
to the usual QED bremsstrahlung and follows from the coherence of the soft 
gluon radiation and the angular ordering condition which limits the angles 
of the secondary particles in (|17D by the (typically small) emission angle -t^i 
of the first gluon. 

If the running of the coupling is taken into account, the computation 
already of the leading contributions yields quickly rather lengthy formulae 
with increasing order q. For the first two moments we find 

niX^, Y) = P^iY + A) In - X^}, (21) 

c^'\X^, Y) = (3'{[X{AY + 3A) - X^{X^ - 2Y)] In 

- (3A + AY - 3X^/2) } (22) 

^ 2(3^Y{{X^ + 2A) In ^^1±1 _ 2X^}, (23) 

A 

where the latter approximation applies for ^ Y . Our essential results 
about the behaviour of correlations for X^ — * coincide with the results for 
fixed coupling as expected by the general argument. 

The different approximations for the second moment F*^^^ discussed here 
are compared in Fig. ^ One can see that the running results yield lower 
values for the moments with increasing X^ but all results approach the same 
linear behaviour with the same slope for X^ — > 0. 

It is also interesting to note that the results ( p!^ - |20| ) for the normalized 
moments are independent of the coupling 7q and therefore also independent 
of the colour of the primary parton, i.e. the normalized moments in quark 
and gluon jets approach the same limiting behaviour. The reason is that in 
this limit all gluons are emitted from the primary parton; this is different 
from the case of global moments, where only the first gluon in the cascade is 
emitted from the primary parton within the DLA. 
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3.2 Spherically cut phase space 

Next we consider the multiplicity distribution in the reduced phase space 
with cut in the modulus of momentum k < fc'^"*, but still limited by the jet 
opening angle 9. 

In the present approximations the initial momentum P enters the final 
results only through the boundary A; < P if the leading particle is disregarded. 
This implies that we can easily construct from the previous subsection the 
results for the combined cuts k± < fc^"* and k < fc'^"* by replacing P by A;'^"* 
or y by X = ln(A;'^'"*0/Qo) in the equations for moments such as (|20D. We 
note that the largest deviations from the DLA are expected to come from 
the large momenta k where the energy-momentum constraints and also the 
explicit form of the splitting functions become important. Therefore, a DLA 
formula like ( [20| ) becomes more realistic for restricted momenta k < A;^"*, i.e., 
if we substitute X for Y keeping X smaller than the global scale Y at the 
same time. 

Now the results for the spherical phase space with k < /c^'"* (without 
additional bound on k±) are obtained from the formulae of the previous 
subsection after letting /cj** A;'^"* or Xj_ X and also replacing Y by X. 
For fixed coupling in lowest order we obtain from (plSf) 



The normalized moments K^'^^ are X-independent. This shows the differ- 
ent behaviour of the cut moments in cylindrical and spherical phase space: 
they decrease to zero or stay constant, respectively, if the cut-off is lowered 
towards the kinematic limit. Consequently, as announced in the Introduc- 
tion, soft gluons with limited momentum k have essentially a non-Poissonian 
multiplicity distribution, while those with limited k± are indeed produced 
independently. 



n(X,y)=7o'XV2, 
ir(^)(X,F) = 2''-V,-i/(2g-l). 



(24) 
(25) 
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3.3 Role of angular ordering 

In order to demonstrate the importance of angular ordering for these pre- 
dictions, we arbitrarily drop this condition in the corresponding integrals. 
Let's consider first the integral as in ( ]T7| ) in the case of the cylindrical phase 
space. Taking the second moment as an example, we then obtain in the 
approximation Xj_ <^ Y and a<j fixed 



jQo/P Vi JQo/Si ki JQo/ki V2 ^2 



^ lo^iy (26) 

where the upper limit of the 'i92-integral is shifted from ■i^i to 0. The result 
is symmetric in X± and Y as is the corresponding result for the multiplicity 
([T6|) which applies also to the present discussion; then the normalized moment 
becomes constant 

^ 1. (27) 

On the other hand, if angular ordering is taken into account, one obtains the 
asymmetric result ([T9|) - corresponding to Xj_Y with the stronger X± 

dependence - and finally the Poisson distribution in the limit Xj_ —>■ 0. 

Of course, angular ordering has also consequences for the results in the 
spherical phase space and in general reduces the correlations, for example, 
/^(2) — 1 from (25) as compared to (^7]), however its consequences are less 
dramatic in this case. For the k± cut the primary emission angle is restricted 
towards small angles and then the same holds for the subsequent emission 
angles; in the case of spherical cut the primary emission angle is only limited 
by its maximum, the jet opening angle O, and therefore the restrictions 
for the subsequent emissions are less severe. So the kj_ cut results are more 
sensitive to the constraints from coherence and angular ordering and therefore 
yield the very different results. 
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3.4 Alternative phase space boundaries 

Finally in this section, we want to discuss a different upper bound for the 
momenta which connects to other results. In case of cylindrical cut we have 
considered so far the trapezoidal boundaries Qo < < fc^"* and /i;_l/0 < 
k < P. In this case X± < Y where X± = Y corresponds to the triangular 
boundary of the spherical cut. For comparison we consider now a phase space 
boundary which is given by a second cone with opening angle G' < 0. In 
this case the upper limit of k, at fixed k±, becomes k < k±/Q'. Then the 
angular integrals factorize and, instead of (|T6]) and (p!8| ) one finds simply 

n{X^,AY)=^^oX^AY (28) 
K(^) = lIszL ( (29) 



2q - 1 VAF. 

where AY = ln(0/0') denotes the rapidity difference of the angular bound- 
aries. These results again clearly demonstrate the different roles of longitudi- 
nal and transverse directions: for X± ^ AY the correlations vanish because 
of angular ordering as discussed above and the formulae (28) and (29) cor- 
respond to the results (|16D and ([T9|) for small X± with Y replaced by AY; 
on the other hand, for Xj_ ^ AY the collection of the many produced par- 
tons yields large correlations. In fact, in this configuration the higher orders 
become important with increasing X± and the fractal structure of the par- 
ton cascade yields the power behaviour of the correlations ( "Intermittency" ) 
P1| , P^ , P^ . Eq. (29) nicely demonstrates the different trends K^''^ and 
K^''^ CO for X± and AY respectively. 



4 Multiplicity moments in higher orders 

In this section we derive the integral evolution equations for factorial mo- 
ments in the cut phase space. These equations will then be used to generate 
in the systematic way higher order corrections to the required accuracy. The 
complete solution obtained in this way gives cut moments for the full range 
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of the cut-offs < X, X± < Y. It reduces to the Bom approximation at the 
threshold, thus justifying simphfications used in the previous Section. For 
simphcity, only the gluon jets will be considered in this section, hence we set 

Ca = 1- 

All evolution equations, which we seek for, follow from Eqs.(7,8) together 
with ( pUD which is also valid in the restricted phase space. To begin we first 
simplify Eqs.(7) employing the pole approximation This essential step 
consists of saturating the angular integration in (7) over the intermediate 
parent momentum K by the leading singularities of the integrand. This 
amounts to approximating K \\ ki in all nonsingular terms and integrating 
only over the smallest angle Qxki, ki = 1 . . .n. This leads to the following 
equation for the single parton density differential in momentum k and angle 



^ 14 



k^ §JkKJ3si Qj 

(30) 

For the fully differential two-parton connected correlation function 

rg^(^i, ^2, h, k^) = )(^9i, ^,2, h, k,, p) + T^?{d2, ^,2, k2, k,, p), 

(31) 

we obtain 

rf\^i, ^12, ki, k2, P) = Mp{k,)p'^'\k2, ^2, K) + ^ r '^iliK'd, 

Vi Jk> K 

\ dQKkJ'\QKk,,k,,K)p^'\^,2,k2,K) + 

r dQKk,Tf\QKk,,^i2, K h, K)\ . (32) 

J^12 J 

with the polar angle of particle i, the relative polar angle ■})i2 of two par- 
ticles and fc> = maxjfci, ^2}- Note the different bounds on Oxki in the two 
terms. They follow from the different singularity structure of the product 

term and of the connected correlation function. The splitting (^) is a con- 

(2) 

sequence of the pole approximation - results from the saturation of the 
d^K integration in Eq.(7) by the K \\ ki configuration. 
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Analogous simplifications of the evolution equations can be made for 
higher order densities. 

4.1 Evolution equations for momentum cut moments 

Integrating Eq.(|30D over the spherically cut phase space of the single parton 
in a jet (P, B) gives for the average multiplicity 

(33) 

The structure of this equation can be simply understood from the general 
rules of the angular ordering and transverse momentum limitations. Since 
the momentum of a child parton is limited by /c'^"*, their emission angles are 
bounded from below by Q^/k'^^*. Because of the angular ordering, this also 
limits from below the emission angle d of the intermediate parent K. The 
remaining phase space of {K^ d) can be conveniently split into K < A;^"* and 
K > fc'^"*, c.f. regions (I) and (II) in Fig. ^. Parents in the first region 
contribute to the global multiplicity at the scale k'^'^^ and parents from the 
second region give rise to the second term of Eq. (|33D. In Appendix A we 
derive an analogous equation for the second cumulant, by following the above 
steps in detail. 

Similarly, one derives the evolution equations for higher moments. We ob- 
tain for the cumulants c^'^^k^'^^, PQ) = c^'^\X, Y) in the logarithmic variables 
X = In (fc^"*e/Qo) and Y = \n (PQ/Qo) 

Jo 

rX rY~X+x 

+ dx dy^l{y)f^''\x,y). (34) 

JQ Jx 

This, together with Eqs.(|10|,|lT]) , uniquely determines all multiplicity mo- 
ments in the spherically cut phase space. The lowest order contributions are 
given recursively, in g, by the first integral with the /^''"^^ term only. Per- 
turbative expansion in the order 0{^q^) can now be generated recursively 
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in q and k by integrating the 0(7g ) expansion according to Eg. (|3^ . For 
constant as this was done for the first four moments up to terms 0(7o^) with 
the results shown in Fig. ^. A sample of explicit expressions is given in the 
Appendix B. 

In the global limit (X Y) the second integral disappears, and the 
above equations reduce to the known equations ( ]TTD for the global moments. 

Since the lowest order is also given by the first integral, the Born ap- 
proximation for the normalized moments is independent of the maximum 
virtuality PG and equals to the global moments at the scale k^'^^Q. Conse- 
quently, the normalized moments are constant in this approximation and are 
determined by the coefficients (|T3p. In the higher orders they aquire a mild 
X dependence due to the additional radiation from parents which are faster 
than the cut-off fc'^"*. 

As already mentioned, while the higher order corrections are of course 
important for the unnormalized moments, they largely cancel in the ratios 
As a consequence, the Born approximation discussed in the previous 
chapter reproduces rather reliably the properties of these moments, especially 
the difference between the spherically and cylindrically cut moments. 

This scheme applies for the running case as well, and higher order 
corrections can be generated by the same, albeit more tedious steps. 

4.2 Cylindrically cut moments 
4.2.1 Multiplicity for constant 

Integrating Eq.(|30|) over the parton momentum k at fixed k±_ gives the evo- 
lution equation for the inclusive fc^ distribution 

^"^ rP rIK rk± K 
p{k^, PQ) = bik^, PQ) + ^ ^ dK^p{K^, KQ^) 

k± Jkj^/e K JQo k± 

(35) 

with the Born term 6 = 7q In {PQ/k±)/k±. Note that the kinematics of this 
equation is very different from that for the absolute momentum distribution. 
While in the latter the rotationally invariant momentum of a child parton is 
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the same for the P-jet and for the internal K-jet, in the former the transverse 
momenta k± and k_l refer to different axes, see Fig. This leads to the more 
complicated structure of Eg . (|35|) . In particular, the virtuality of the internal 
K-jet depends now on k±. In fact the effective angle = B|-|- should be 
interpreted as the angle between the momentum of the parent K and the 
parton emitted at the maximal angle B which is allowed. In this variable, 
Eq.(^) has a simple interpretation. 



Integrating Eg . (|35|) over k± up to the transverse cut-off A;^ , we obtain 



for the average multiplicity of partons with limited kj 



n(A;l"*,Pe) = nb(A;l"*,P0)- 



Qo/e K Jqo/k -d 

rP dK /•© d-d i9 
ll -^n{kr%m. (36) 

where fik denotes the Born contribution, Eg. ([l6|) . Similarly to the spherical 
case, the phase space of the intermediate parent is divided into two regions. 
In the first integral the maximum virtuality of the K-jet is limited by fc^"*, 
therefore all partons in that jet are below cut-off, hence the global multiplicity 
contributes. In the second integral some partons may exceed the external cut- 
off fc]"*, thus only cut multiplicity contributes. The cut k™"^ = (k'^^/Q)^ 
corresponds to the hardest parton compatible with the cut-off fc^* emitted 
from the parent K at angle i). The above eguation reproduces the exact 
solution which can also be obtained, in the case of constant a^, by direct 
integration of the known fully differential distribution. 

4.2.2 Relation between spherically and cylindrically cut moments 

We have seen that the natural variable appearing in the k± cut case is the 
angle between the hardest parton in an original (P, G) jet and the interme- 
diate parent K. One can however reorganize the calculation such that the 
phase space of the intermediate parton K is parametrized by its momentum 
and the angle Qpx at which it was emitted from the original parent P. This 
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formulation provides very useful relations between both families of cut mo- 
ments, valid also for running a^. It will be used to generate cylindrically cut 
moments to higher orders. 

We begin again with the first moment, i.e. the average multiplicity writ- 
ten as, c.f. Fig. ^ 

nik"^', PB) = nik^f) + dkp^^\^,k,P). (37) 

Due to the independence of the fully differential distribution on the opening 
angle 0, the integration over the small angles of the emitted partons Qo/P < 
^ < k'J^^/P results in the global multiplicity at the scale fc^"* (region A in 
Fig. 1^). The second term represents contribution from the region B. We 
use now the Eq.(^) to express the last term by the contributions from the 
intermediate parents. Rearranging the integrals gives after some algebra 

n,yiiX^,Y) = n{X^) + {Y - X^) f''\',iy)[l + n{y)]dy 

Jo 

+ r (Y- yh'oiy)^sph{Xu y)dy. (38) 

where the subscripts distinguish different cuts, and logarithmic variables are 
used. Here and in the following the quantities without subscripts and one 
argument refer to the global quantities. This equation allows to calculate 
fc^-cut multiplicity if the momentum cut result is known. Using the constant 
Us solution of Eq.(|5BD we have reproduced with the aid of (^Sp perturbative 
expansion of ficyi to 0(70"*^^). 

This approach generalizes for the moments of arbitrary order. We begin 
with the q = 2 case and display explicitly the relevant two-parton phase 
space for the global cumulant 

s^\pQ)=2r dA r r dk,r dk2T^^\A,A2,k,,k2,p), 

JQo/P JQo/P JQo/^12 

(39) 

where F^^^ defined as in Eq . (|3l|) . The factor 2 comes from the symmetric (1 ^ 
2) configuration. Fully differential connected correlation functions in the 
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leading logarithmic accuracy depend only on one angle relative to the original 
parent momentum, e.g. Direction of the second particle enters only via 
the relative angle i)i2. Moreover for ■}}i2 > i^i the connected correlation 
function vanishes - yet another manifestation of the angular ordering . 



Finally, the fully differential connected correlation function is independent of 
the global opening of the original jet as in the single parton case. The last 
property allows us to split the cut cumulant into two parts as previously 



/ dkj dk2Tf\A,A2,h,k2,P). (40) 

Since the leading singularities are generated from the region where in fact 
^12 ^ the measured transverse momentum of the second parton k2± ~ 
k2^i and this gives the upper limit of the ^2 integration in pO|). Inserting 
now Eq.(|32D under the integral in pO| ) and rearranging various integrals gives 
the final expression for the cylindrically cut cumulants in terms of the global 
and the spherically cut factorial moments. Details of this calculations are 
presented in Appendix C. 

Similar steps allow to derive analogous relations for cumulants of arbitrary 
order q. We obtain 

{Y - X^) il{y)[qf^'^-'\y) + f^'^\y)]dy + 

r {y-yH{y)fit{x^.y)dy. (4i) 

In fact Eg. (p8|) is the special case of this relation for g = 1 since f^^^ = 1. 

Together with Eqs. (p^ this equation uniquely determines all moments 
with limited k±. We have used these equations to generate algebraically the 
perturbative expansion of the first four moments up to the eight order in 7q 
for constant a^, see Fig.^ and Table 4. We quote only the first nontrivial 
example - the 0{a^) correction for c'-^-* 

4S,3(^±, y) = ^^1 - \xiY + Ixiv. (42) 
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This result was confirmed independently by the direct integration of the 
resolvent representation for p^^^ [l^ and also from the first explicit iteration 
of Eq.(7). 



4.3 Properties of the higher order results 

Fig. 1^ shows spherically and cylindrically cut moments up to 0(af) in the 
full range of the cut-off. First three normalized moments are displayed. 
All properties discussed earlier are fully confirmed by this calcuation. For 
small cutoffs all cylindrically cut moments tend to unity while those with 
limited momentum have different limits corresponding to the non-Poissonian 
multiplicity distribution of soft gluons. The perturbative expansion is rather 
poorly convergent for the unnormalized moments which, being polynomials 
in the cut-off and Y, grow faster for higher order k. However this growth 
is largely canceled in the ratios. As a consequence, perturbative expansion 
is remarkably stable for the normalized moments F^'^\ This is clearly seen 
in Fig. 1^ where different lines, corresponding to different maximum order 
kmax included, group together for each q and each type of the cut-off. In the 
global limit X_i_,X = Y spherically and cylindrically cut moments are equal. 
However this common value, which is in fact given by the global moment, 
does vary with the order of the perturbative approximation. This is seen 
in Fig. ^ as well. This dependence is very weak and corresponds, for large 
Y, to the difference between the coefficients which determine the threshold 
behaviour (see Eq.(0)) and the asymptotic behaviour (see [l^), respectively, 
of the global moments. 



5 Monte Carlo results 

We have studied the behaviour of the cut moments within a Monte Carlo 
program in order to check our analytical calculations with a more complete 
numerical method, which is important in the test of the LPHD picture. In 
view of the similarity of its scheme with our analjd^ical framework, we have 
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chosen the ARIADNE program|T9[, in which the perturbative phase is ter- 
minated by a cutoff in transverse momentum and, in its defauh version, 
followed by string fragmentation. In deviation from the default value, the 
cutoff Qq for partons can be chosen close to the QCD scale A as required in 
the LPHD picture, which is not easily possible in other popular QCD Monte 
Carlo models. 

For our applications we have then switched off the hadronization phase 
within ARIADNE and directly compared the result at the end of the per- 
turbative evolution with our theoretical predictions. We have reset the 
values of the parameters Qo and A in ARIADNE in order to directly re- 
produce the energy dependence of the average multiplicity of all (charged 
plus neutral) hadrons according to the LPHD picture. We estimated the 
full (charged plus neutral) multiplicity as 3/2 the measured multiplicity for 
charged hadrons and we fixed the overall normalization factor which relates 
partons to all hadrons to Kaii = 1, as suggested by previous studies on the 
average multiplicity pO|0. Table |l] shows a sample of parton multiplicities 
obtained by ARIADNE for different parameters Qo and A. By comparing 
the parton results with the experimental data on all hadrons, one concludes 
that the best set of parameters is given by Qo = 0.2 GeV and A = 0.015. It is 
remarkable to notice that the same value of A has been obtained in previous 
studies on jet and particle multiplicities pOU , whereas the present value of Qo 
is slightly smaller than in the previous case. This difference can be related 
to the different approximation schemes applied. 

As a further test of LPHD, let us compare the predictions of ARIADNE 
using the new parameters with experimental data on factorial moments in 
one hemisphere as measured by the OPAL Collaboration [^. This is done in 
Table |^, which clearly gives further evidence in support of LPHD using the 
Monte Carlo results at the parton level. In this comparison we have assumed 



^There is a correlation among the parameters Qo and Kaii; the chosen values provide 
the best description of the average multiplicity not only at LEP-1 cms energy but also at 
lower cms energies 
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v/i = 91.2 GeV 


^/s = U GeV 






=31.0 ±1.1 


flail = 13.95 ± 0.65 




A 


n„ 

y 


n„ 


0.6 (def.) 


1 (def.) 


8.21 


4.2 


0.27 


0.1 


19.2 


9.7 


0.27 


0.01 


24.8 


12.5 


0.23 


0.01 


29.4 


13.4 


0.2 


0.015 


30.2 


13.8 



Table 1: Predictions by ARIADNE at parton level for the total average 
multiplicity fip in e~^e~ annihilation at ^/s = 14 and 91.2 GeV for different 
values of the parameters Qo and A, compared with the experimental results 
for all hadrons fiaii (taken as 3/2 the value for charged hadrons). 

that the normalized higher order global moments, Fq, are equal for charged 
hadrons and for all hadrons. 

After having fixed the set of parameters of ARIADNE which are consis- 
tent with the picture of LPHD, let us now study the predictions for the cut 
moments at parton level; the dependence of the A;^"* ^ind /c™* moments on 
the cut parameter in one hemisphere defined via the thrust axis is shown 
in Fig. ^. According to our analytical calculations, we expect that the k± 
moments approach unity as In k±/Qo for k± ^ Qq. The cut moments pre- 
dicted by ARIADNE show indeed such behaviour for small values of the cut, 
^cut ^ ^ GeV; however, at very small A;]"* —>■ Qq, the moments do not reach 
the predicted Poissonian value of 1, but saturate at a value larger than 1. 
Such effects are expected from threshold effects close to the infrared cutoff, 
in particular, from the boost needed to relate the overall rest frame of the 
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n 


ir(2) 


ir(3) 


ir{4) 


Partons 


30.18 


1.09 


1.30 


1.68 


Hadrons 


30.3 ±0.6 


1.0820 ±0.0052 


1.275 ±0.019 


1.637 ±0.047 



Table 2: The mean multiplicity n [up or naii) and the first three normalized 
factorial moments in one hemisphere defined through the thrust axis pre- 
dicted by ARIADNE at parton level with parameters Qq = 0.2 GeV and A = 



0.015 and experimentally measured by the OPAL Collaboration |^T|. Errors 
on data are the sum of statistical and systematic errors. 



collision to the single dipole rest frame where each new emission takes place. 

The k'^^^ moments in Fig. rise to rather large values at intermediate 
scales and bend at small A;'^"*; the important point here is that they reach 
finite values much above unity, showing that in spherically symmetric phase 
space there is no Poissonian regime at very small momentum cuts. Both types 
of moments approach the same values for large /c^"* ~ fc'^"* corresponding to 
the global moments as it should be. 

To investigate in more detail the behaviour of the k± moments in the 
kinematic region close to Qo, let us now examine a new simulation with the 
same value of A but the larger value Qo = 0.4 GeV. Results of this run are 
shown in Fig. |^ for both A;]"* and /c™* moments. In this case, in which the 
cascade is stopped earlier, the situation is clearer: the k± moments approach 
the Poissonian value one for fc^"* ^ Qo theoretically expected, and stay 
there below the cutoff. The /c^"* moments saturate, on the other hand, at 
values above one as before. 

From these results we conclude that the more complete Monte Carlo 
calculation confirms the essential features of our analytical calculation in the 
region of small cut variables: for the kf'^-cut the moments vanish linearly as 
ln(A;]"*/(5o), whereas for the k'^'^^-cut the moments approach a constant value. 



21 



However, the slope and the constant values are different from the analytical 
predictions and also the dependence for large values of the cut momentum. 
We attribute these differences to our simplified DLA calculations which does 
not take into account energy conservation constraints and the nonleading 
parts of the parton splitting functions. 

Finally, we show in Fig. ^ the predictions for the cut moments both with 
cylindrical and spherical cuts obtained at hadron level by running ARIADNE 
with default parameters for e~^e~ annihilations at cms energy of LEP-1. Data 
refer again to particles in one hemisphere, defined through the thrust-axis. 
Electrons and positrons have been subtracted to avoid the lepton contamina- 
tions of the hadronic signal; by keeping the electron pairs, a peak at very low 
values of the cut momentum would appear, masking the eventual physical 
signal. This result is consistent with an old observation] 2^ done within the 
study of intermittency phenomena that e^e~ pairs from vr'^ decays strongly 
affect the experimental results in very short-range correlations. 

As expected, for large values of the cut, the two moments go to a com- 
mon value corresponding to the global moments of the multiplicity distribu- 
tion without any cut; the obtained results are close but up to about 10% 
smaller than the experimental values for the global moments directly mea- 
sured by the OPAL Collaboration [^. For small values of the cut, the two 
moments show a different behavior; the cylindrical symmetric cut moments 
for A;^"* ^ 5 GeV first tend towards unity, in agreement with our theoretical 
expectations; however, they rise again for very small values of the cut. This 
effect does not occur at the parton level and therefore has to be associated 
with hadronization corrections, in particular production and decay of reso- 
nances|. The spherical symmetric cut moments show a step around a few 
GeV's which separates two plateaus. The depletion at large fc"^"* could be 
interpreted again as the result of energy-momentum constraints not included 



^In order to check this point, we have swichtcd off the perturbative cascade within 
ARIADNE and found that all moments approach unity for fc™* > 1 GeV, while the peaks 
at small values of fc™* remain. If resonances are not allowed to decay, all moments become 
more reduced at small fc™* , but small peaks are still visible. 
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in our DLA calculations. It is remarkable to notice that the moments satu- 
rate and do not fall down for small values of the cut momentum, very much 
like in our analytical calculation. However, in view of the parton level results 
for the same quantity we consider this coincidence as accidental. 

By comparing the partonic and the hadronic predictions, one sees that the 
LPHD is approximately confirmed within ARIADNE at the level of global 
moments and for large values of the cut momenta (right hand side of the 
plots) but is clearly violated for soft particles (left hand side of the plots). 
So we conclude that the parton level predictions show indeed the peculiar 
features predicted analytically within the perturbative approach, namely 
the near Poissonian production of soft particles as consequence of coher- 
ent multigluon emission, whereas the same features are softened and partly 
washed out in the model with string fragmentation. It will be interesting to 
find out the experimental situation. 

6 Conclusions 

We suggest analysing multiparticle production in restricted phase space re- 
gions with variable cut A;^* in the transverse momenta of the particles, and, 
for comparison, with variable momentum cut A;™*. We have derived the 
evolution equations for the multiplicity moments within the DLA of pertur- 
bative QCD, from which their dependence on the jet virtuality and on both 
cut variables can be obtained. Results for moments up to order q — A are 
given explicitly in case of fixed coupling by the first terms of the perturbative 
expansion. 

The normalized factorial moments are found to approach a linear be- 
haviour for small values of the transverse momentum cut with F^'^^ 1 for 
fcj"* — > Qo corresponding to a Poisson distribution, as the parton showering 
is uniformly suppressed at small k± for any rapidity. On the other hand, for 
the momentum cut, this is not the case and the moments approach a finite 
value for the minimal value of A;'^"*. 
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The DLA takes into account only the leading singularities but respects 
the soft gluon coherence which is implemented through the angular ordering 
prescription. We have shown that this property is responsible for obtaining 
a Poisson distribution in the soft limit. As a check of our analytical results, 
we have compared them with the results of a parton Monte Carlo program 
which fully takes into account the constraints from energy-momentum con- 
servation and the complete parton splitting functions. The factorial moments 
for decreasing A;™* continuously decrease to values close to unity, contrary to 
the case of the momentum cut, as expected from the analytical calculations. 
At the quantitative level, however, the DLA results show considerable devi- 
ations from those of the more realistic Monte Carlo calculations, as already 
observed in the case of global moments. 

It will be interesting to find out to what extent the experimental data 
follow the perturbative results as one might expect from previous successes 
of the LPHD picture. Final state interactions, in particular resonance pro- 
duction, could severely disturb the perturbative results, especially in the soft 
region, as it is suggested indeed by the results of a Monte Carlo program 
which implements a model of the hadronization phase. By the proposed 
studies the limitations of the perturbative predictions could be explored in 
genuine multiparticle correlations. 
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Appendix A 



We derive here the evolution equation for the second cumulant in the spher- 
ically cut phase space following in detail the steps outlined in Section 4.1 . 
Integrating Eqs.(7,8) for g = 2 gives 

+ / d^KMp{K) [ d^ki [ d^k2p^Piki,k2), (43) 

JTi+il+ui Jki<k'=^t Jfc2<fe'="t 

where the order of integrals in the second term was interchanged. The full 
phase space of the intermediate parent K is displayed in Fig.|^. As in the 
one-parton case, only regions I and II contribute since both momenta of final 
partons are limited from above by the cut-off k'^'^^. This in turn restricts from 
below their emission angles and consequently, due to the angular ordering, 
the emission angle of the parent K. We then split explicitly the remaining 
integral into regions F/ and F//. Simplifying also the first term gives 

c^'\k^-\P,Q) =2 I ^ f ^7o'(Mi)r^(A:i^i) 

rB di9 rP dK 
+ 1 -J ^^',{m)f^\k^'^\K,^). (44) 

In the first term, the k2 integral gave just the global average multiplicity at 
the scale ki-t} since the restriction k2 < k^^^ is weaker than k2 < ki < k^"^^ 
which is satisfied for this subprocess anyway. Similarly in the second term, 
momentum of the parent K in the region F/ is smaller than the cut-off, 
hence the phase space of the final partons is not restrited by A;'^"* but by the 
momentum of the parent K < A;'^"*. Therefore the global moment at the 
scale Kd results. On the other hand in the third term a parent K is harder 
than the cut-off and phase space of the final partons is indeed restricted 
by A;'^"*. Consequently integrating over ki and k2 gives the cut moment. 
It is natural to introduce the logarithmic variables Y = ln(PB/Qo) ^ = 
ln(A;'="*0/go)- Then c(2)(A;'="*, P, 6) = c^'^\X,Y). Finally we observe that. 
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due to the evolution equation for the global moments, Eq.(^T]), the first two 
terms combine into the global cumlant at the scale k'^'^^Q, hence 

c^'\X,Y) = c^'\X) + f^dx dy^l{y)f'\x,y). (45) 

Jo Jx 

Appendix B 

In this Appendix we collect a sample of higher order expressions for the cut 
moments. The perturbative expansion for a generic (spherically or cylindri- 
cally cut) moment reads 

flLiZ,Y)=Y.lo f!f\z,Y), (46) 

k=q 

where the cut-off Z = X for moments in a spherical phase space and Z = 
Xj_ for a cylindrical cut. Our results for the coefficients flp\^j^{X,Y) and 
fcyi,k{'^±^^) shown in Table 3 and Table 4 respectively. The formulas 
for the spherically cut moments were generated (by Mathematica) recursively 
in k and in q from the evolution equation, Eq.(|5^) coupled with Eq.(|TU]). 
The formulas for cylindrically cut moments were obtained by integrating 
spherically cut moments according to Eq. (^T)) . As one of the consistency 
tests we have checked that the results obtained from Tables 3 and 4 for the 
normalized moments F^'^'^ = f^f' /(fi^^ Y, in the leading order of X and 
agree with those from Eqs. (25) and ( [2D| ) respectively. 
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k 


fsph,ki-^-> ^) 




1 


2 


1 


2 


8 6 




3 


x^Y 1 x-^y^ 

72 30 ' 48 




2 


3 


2 


3 


X6 7X^Y 
18 30 




4 


240 ~'~ 630 ' 720 




3 


17X6 
60 


3 


4 


79XS 1 43X''y 
1120 ' 140 




5 


23X^0 1 647X''y 1 323X«y2 
2800 ' 10080 ' 2240 






4 


31X« 




105 


4 


5 


737X10 1 412X9y 
7560 ' 945 




6 


76603X12 4757X"y , 6161Xi''y2 




4989600 ' 37800 ' 21600 



Table 3: Perturbative expansion of the spherically cut moments 
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k 






1 




1 


2 


8 3 ' 4 




3 


Xl x'lY x|y2 x^ys 

72 20 16 ' 36 




2 


^^^^ + XlY^ 


2 


3 


xl xlY 3x4y2 x3y3 

36 3 4 ' 2 




4 


x» iix^y 89XKy2 loix^y:^ I7x4y4 

90 120 ' 360 360 ' 144 




3 


xe 3X5 y ^45.2 
12 10 2 T^.L-' 


3 


4 


29X8 83X7 y 11x6 y2 5 3 3X4 y4 

480 ' 420 ' 60 ^ -L ' 4 




5 


551X]0 589X9y 97X8y2 ii7xXy3 857X6y4 I3x5y5 
33600 10080 2240 ' 280 1440 ' 48 




4 


Xl , 2X1 y 13X6 y 2 , ^4 

15 "^21 15 + ^-L^ 


4 


5 


iix]o 383X9 y 169x8 y2 iix^y^ 5x6 y* ^5^^ 

280 3780 ' 210 15 6 + 




6 


1087X12 14941Xiiy 683X]"y2 27017X^y3 187X8y4 37XXy5 35X6y6 
151200 ' 277200 1800 ' 45360 ' 2520 45 ' 72 



Table 4: Perturbative expansion of the cylindrically cut moments 
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Appendix C 



We will analyse in detail the steps leading to the evolution equations for the 
cylindrically cut moments, Eqs. (|38| , ^TD . The cut multiplicity is defined, c.f. 
Fig. I 

= / dkp{^,k,P)+ d^ dkpi^,k,P). {A7) 

JQo/P JQoM Jk'j^yp JQo/i} 

The first term (region A) is just the global multiplicity at the scale /c^*. In 



the second term we use Eq. (|30|) . 

r© d^ r^T'/^ 



dk 

1^ 



n(A;r,P,e) = n(A:r)+ / ^ / 

Jk^^^yp if JQo/ 

\^-^+r^ll{mf dQKkPiQKk.KK)]. (48) 

{ k Jk K JQo/k J 

which after changing the orders of dk and dK integrals gives 

nikT\ P, e) = nik^) + f ^ / " ^l'om + 

Jk'j^yp V JQo/^ k 

— / — 7o(^^) / dk / dQKkPiQKk, k, K) + 

di9 rP dK rk'^^y-d 

— -rrloim dk dQKkPiQKk,k,K). (49) 

kcut/p ^ Jfc^*/^ K JQo/-d JQo/k 

In the third term integrations over k and Qxk extend over the full phase 
space of the K jet, hence they give global multiplicity. In the fourth term 
the child momentum is restricted, therefore the momentum cut multiplicity 
results 

n{kr,p,e) = nikr) + ^ ] iriim+ 

Jk^^^P V JQo/S k 
fkr/^dK 2,,^ ON /r^oN 

kcnt^p ^ jQo/-d K 

/V,? rP (]K 

^ -jTllimnsAkTM.K,^) (50) 
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which when transformed to the logarithmic variables gives Eg. (|3^) . 

Let us now study the second cumulant, Eqs. (|39D . Inserting Eg . (1321) in 



(]40|) gives 

.(2) 



2 / , l-7o(^i^i) 

d^i2 r dk2p^^\k2,^i2,ki) + 

2 / -r^ / —jn(K^i) / d^i2 
Jk^yp ifi Jqo/-&i K Jqo/k 

r dki deKkJ'\&Kk,,ki,K) 

r dk2P^'\A2,k2,K) + 

© d§l fP dK 2/..n N f^^ 



2 7o(^^i) / 

dki dQxk^P^ '{QKki.ki, K) 

'o/^i2 JQo/ki 

dk2p^'\^u,k2,K)+ 

Qo/^12 



'12 



j-e f^TI^^ dK ^^^^^ ^ f^i ,^ 



dk, I dk2 I ' deKk,T^^\<dKk,,A2,k,,k2,K) + 



,© d^i [P dK 2..^o ^ f^^ 



cut 



'12 



dk, dk2 dQKk,T^^\QKk,,^12,k,,k2,K). (51) 

?0/'fl2 JQo/i}i2 J^12 

In the second term the integrals over i)i2 and fcihave been interchanged. Since 
this term describes the process P k, ^ k2, a. momentum k2 is limited in 
fact by k, < kf'^/i}, and not by the cut-off k^^^ /d, itself. In the remaining 
four terms the K integrations have been pulled in front of the k,, ^2 and 'di2 
integrations. It is then convenient to split the full range of the intermediate 
parent momentum K according to the cut-off kf-^/^,, similarly to the average 
multiplicity (g=l) case discussed above. For K below the cut-off the inner 
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integrals cover the full phase space of the K jet, while for K above the cut- 
off only the spherically restricted integrals over children momenta ki and ^2 
occur. All these inner integrals can be cast into the global or spherically 
restricted moments. We have 

cSi(A^r,p,e) = c(2)(A:r)+ 

r® di9, rP dK 

/■e d§i fkr/^^ dK 2,^^, , (2^^r.a^ 

^l'oiKA)c%{kr/^uK,A)- (52) 

Finally, introducing logarithmic variables, and after combining cumulants 
and squares of the average multiplicity into factorial moments, we obtain the 
g = 2 counterpart of Eq. (|l|) . 
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Figure Captions 

Fig. |l|. The normalized factorial moment F'-^-' = 1 + K'^'^^ for different 
approximations in leading order. The lowest curve represents the result (22) 
for running a^, the next one the small X±_ approximation (23). The upper 
two curves represent the fixed as result following from (|TB|) and the linear 
approximation (pO|). In all cases F = 5.1 and A = 0.015. 
Fig. ^ Phase space of an intermediate parent for the momentum cut mo- 
ments, c.f. Eqs.((|33| ) ,(|43| ) ). Parents in region I generate full jet with virtuality 
f^cutQ^ the ones in the region II contribute to the cut moments, while region 
III does not contribute due the angular ordering. 

Fig. |3|. DLA predictions for the cut-off dependence of the first three normal- 
ized moments F^'^J^^^, q = 2,3,4 for Y = 5.7. Both families, i.e. spherically 
(sph) and cyllindrically (cyl) cut moments, are shown. They coincide in the 
global limit (X or Xj_ = Y), but have distinctly different threshold behaviour. 
Different lines which describe one moment correspond to different order of 
the perturbation theory included kmax, kmax = <?, g -|- 1, 8. 
Fig. ^ Kinematics relevant for the evolution equations at fixed transverse 
momentum, Eqs.((p5|),(p^)); 'd = 0kj_//c_l is an angle between a parent K 
and a parton emitted at the maximal angle G. 

Fig. ^ Kinematics relevant for the splitting the cylindrically cut moments. 



c.f. Eqs.( (|37|) ,(40)). Partons in the region A form a complete jet with virtu- 
ality fc^"*. 

Fig. a: cut moments of order 2 (diamonds), (crosses) and 4 (squares) 
in one hemisphere defined through the thrust axis as a function of fc^"* as 
predicted by ARIADNE at parton level with parameters Qo = 0.2 GeV and 
A = 0.015 at ^/s = 91.2 GeV. b: same as in a, but as a function of /c'^"*. 

Fig. 0. Same as in Fig. §, but with Qo = 0.4 GeV. 

Fig. Same as in Fig. but at hadron level with default values of the 
parameters and string fragmentation. 
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Figure 3: 
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Figure 4: 
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